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Abstract 

An analysis is given of the effects of non-universal soft SUSY break¬ 
ing masses in the Higgs sector and in the third generation squark sector, 
and it is shown that they are highly coupled. Analytic expressions are 
obtained for their effects on the parameters p, ruA and on the third gen¬ 
eration squark masses. Non-universality effects on dark matter event 
rates in neutralino-nucleus scattering are analysed. It is found that the 
effects are maximal in the range < 65 GeV where the relic density 
is governed by the Z and Higgs poles. In this range the minimum event 
rates can be increased or decreased by factors of 0(10) depending on the 
sign of non-universality. Above this range Landau pole effects arising 
from the heavy top mass tend to suppress the non-universality effects. 

The effect of more precise measurements of cosmological parameters on 
event rates, which is expected to occur in the next round of COBE like 
sattelite experiments, is also investigated. Implications for the analysis 
for dark matter searches are discussed. 

1. Introduction 

Although there is currently a great deal of evidence for the existence of dark 
matter in the universe |l|], the nature of such dark matter is still uncertain. 
COBE and other experiments have provided further insights into the na¬ 
ture of dark matter which strongly suggests more than one component to its 
structure]^. Thus the simplest types of models to fit the power spectrum seen 
by COBE and other experiments require two components, a hot component 
and a cold component. Several other possibilities have also been discussed 
such as those using a cosmological constant. 

In our analysis here we shall adhere to the simple two component picture 
of dark matter, with a hot component (which could be a light neutrino with 
mass in the range of a few electron volts), and a cold component which we 
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assume is the lightest supersymmetric neutral particle, the neutralino. (For 
a review of supersymmetric dark matter see ref.3). For the purpose of our 
analysis we assume a mix of hot and cold dark matter in the ratio of VLhdm '■ 
^CDM = 1:2 (consistent with the COBE data), where hlj = Pi/Pc where pi 
is the matter density in the universe due to matter of type i and pc is the 
critical matter density needed to close the universe. Assuming an inflationary 
scenario with hi = 1, the condition that Qb < 0.1 as implied by the Big Bang 
Nucleosynthesis, and h in the experimental range 0.4 < h < 0.8 where h is the 
Hubble parameter in units of 100 km/s Mpc, one hnds in the range 


0.1 < % < 0.4 (1) 

is the quantity computed theoretically and thus the above condition acts 
as the primary constraint on the nentralino dark matter analysis. We shall 
discnss in the latter part of this paper the effect of constraining the 
range more narrowly. 

The pnrpose of this paper is to analyse the effects of non-nniversal bonnd- 
ary conditions on soft SUSY breaking masses at the GUT scale on event rates 
in nentralino-nncleus scattering, as almost all of the previous analyses, with 
the exception of the work of Ref. (4), are in the framework of snpergravity 
nnihcation using only the nniversal bonndary conditions. 

The ontline of the paper is as follows: in Sec. 2 we give the salient points of 
snpergravity nnihcation which provides the framework of onr analysis. In Sec. 
3 we discuss non-nniversal soft breaking in the Higgs sector and in the third 
generation squark sector and show that their effects at low energy are highly 
conpled. We then obtain solntions in a closed form at the electro-weak scale 
of the relevant low energy parameters that contain the non-nniversalities. In 
Sec. 4 we discuss the constraints that are imposed on the event rate analysis. 
In Sec. 5 we give the basic formulae that enter in the analysis of scattering of 
the nentralino by nnclear targets. In Sec. 6 results of the event rate analysis 
are given and a comparison made between the cases with nniversal and non- 
nniversal soft SUSY breaking bonndary conditions at the GUT scale. In Sec. 
7 we discuss the effects of grand nnihcation on non-nniversal parameters. In 
Sec. 8 we discuss the effects that a more accurate determination of the Hnbble 
parameter, which is expected to occur in the next ronnd of GOBE like sattelite 
experiments, will have on the event rate analysis. Gonclnsions are given in Sec. 
9. New results of this paper are contained in Secs. 3,6,7,8, in Figs. 1-7, in 
Tables 1-5 and in Appendix A. 
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2. Neutralino Dark Matter in Supergravity Grand Unification 

In the analysis of this paper we use the framework of supergravity grand 
unification(SUGRA)0. In this picture one assumes the existence of a su¬ 
pergravity grand unified theory which in its symmetric phase operates in the 
region between the string scale and the grand unification scale. Supersymme¬ 
try is broken in the theory via a hidden sector. In the minimal supergravity 
unification one finds that after breaking of supersymmetry and of the grand 
unified symmetry (where one assumes that the grand unified theory breaks to 
the gauge group SU (3) x SU (2) x f/(l) ) one has a low energy theory (obtained 
by integrating out the heavy modes and the modes of the hidden sector) which 
can be given in terms of just five arbitrary parameters [^“|@]. 

These consist of the universal soft SUSY scalar mass at the GUT scale 
mo, the universal gaugino mass at the GUT scale mi/ 2 , the universal trilinear 
scalar coupling at the GUT scale Aq, Bq which is the co-efficient of the bilinear 
term ^qHiH 2 at the GUT scale, and the Higgs mixing parameter po at the 
GUT scale. In addition there are the GUT parameters Mq and ac which are 
determined by the grand unification condition. 

A remarkable aspect of supergravity grand unification is the breaking of 
the electro-weak symmetry by radiative effects. Radiative breaking is governed 
by the equations § 




2 


tan^(5-l 2^^ “ 2 /i 2 + ^,2 + ^2 


( 2 ) 


where B is the value of the parameter Bq at the electro-weak scale, = 
Tn\j. + Sj where is the running Hi mass at the scale Q ~ Mz and Sj are 
loop corrections ||I0|, pT| . are given by 


=^1 + ml/29{t) (3) 

+ Aomi/2f{t) + mlh{t) - k{t)Al (4) 

Here t=ln(Mg/Q^), and e,f,g,h,k are form factors defined in Ref. (12). On 
using the radiative breaking equations one can determine the parameter p 
from the Z boson mass and one can also eliminate the parameter B in favor of 
tan/? =< H 2 > / < Hi >. Thus in the low energy domain, the theory can be 
described by the parameters 


mo,mg, AtytanjS, sign{^) (5) 

where At is the value of Aq at the electro-weak scale. Further, at the electro- 
weak scale the SU{3) x SU{2) x SU{1) gaugino masses are given by m^ = 
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[piilaG)Ta\i 2 - ( For radiative corrections to this formula see Refs. (13-14)). 
We note that the situation in minimal supergravity is in sharp contrast to 
that for MSSM. In MSSM one has 110 parameters in the low energy domain. 
In contrast, in minimal supergravity unihcation (MSGM) one has only four 
parameters (and one sign) as discussed above. The minimal model has 32 
supersymmetric particles, and all their masses and coupling constants are de¬ 
termined by just the parameters of Eq. (5). 

One of the interesting features that emerges from supergravity unihcation 
is that the model predicts as a consquence of its dynamics which particle is the 
lowest lying supersymmetric particle(LSP). One hnds that over almost all of 
the parameter space that the neutralino is the LSP. Thus with the assumption 
of R parity invariance, one hnds that supergravity unihed theory provides a 
candidate for cold dark matter. The LSP neutralino is an admixture of four 
neutral states, i.e., of gauginos and of Higgsinos Hi,H 2 . Denoting the 

lightest neutralino by Xi we may write 

Xi = riiWs + n2B n^Hi + n4^H2 ( 6 ) 


where the co-efhcients rii are to be determined by diagonalizing the neutralino 
mass matrix given in Appendix A. An interesting property of radiative break¬ 
ing is that it produces the phenomenon of scaling in a signihcant part of the 
parameter space[|l^. In the region of scaling one has >> M|. In this 
region the eigenvalues and the eigenvectors of the neutralino mass matrix can 
be computed in a perturbative fashion by expanding in powers of M^//i 
The analysis shows that in the scaling limit the neutralino is mostly a Bino 
with n 2 > 0.95. Further one hnds that ns, and ni are typically of hrst order, 
i.e., 0{Mz/fj)- These expansions are useful in understanding the relic density 
and the event rates. We add here a note of caution. Inspite of the fact that 
the Higgsino components are generally small their effects on event rates are 
generally very signihcant. Thus large inaccuracies can result from a neglect of 
the Higgsino components. 
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3. Effects of Non-universal Soft SUSY Breaking 

In this section we want to study the ehects of non-universalities in soft SUSY 
breaking on the low energy parameters. While most of the previous analyses 
in supergravity grand unihcation have been carried out with universal soft 
SUSY breaking, the general framework of the theory allows for non-universal 
soft SUSY breaking terms. Thus, for example, under the assumption of a gen¬ 
eral Kahler potential, deviations from universalities can be generated]^. Of 
course the non-universalities are severely constrained, most signihcantly from 
havor changing neutral current (FCNC) constraints. In the analysis here we 
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shall limit the nature of non-universalities to the Higgs sector and to the third 
generation sector. Higgs sector non-universalities are not strongly constrained 
by FCNC and have been discussed in the literature[T3, 1^, ® often in the 


context of certain SO(IO) models to achieve radiative breaking of the electro- 
weak symmetry. In our analysis here we show that the non-universalities in 
the Higgs sector and the non-universalities in the third generation up squark 
sector are highly coupled due to renormalization group effects below the GUT 
scale. It is convenient to parametrise the non-universalities in the Higgs sector 
by 61,62 and the non-universalities in the third generation squark sector by 
^ 3 , 6 i so that at the GUT scale Mg one has 


mH^=ml{l + 6i), mj^^=ml{l + 62 ) (7) 

ml^=ml{l + 63), ml^ = ml{l + 64) ( 8 ) 

where tuq is the universal scalar mass of the hrst two generation masses. We 
shall assume that the remaining sectors of the theory are universal. We shall 
give an analytic solution to the non-universal effects on the mass spectra at 
low energy. We begin by a discussion of the case of the Hi mass and the 
masses of the hrst two generations of sparticles. For we obtain 

mH^=ml{l + 6 i)+my 2 g{t) + ^Sop (9) 

Here the last term arises from the Tr{Ym?) term in the renormalization group 
evolution equations]^. This term vanishes in the universal case because of the 
anomaly cancellation constraint Tr{Y) = 0. However, for the non-universal 
case it is non-vanishing and one has 

rig 

So = Tr{Ym^) = (10) 

2=1 


where all the masses are at the GUT scale and Ug is the number of generations. 
In Eq. (11) p is dehned by 


^ 66^ «i(0)^^ 


( 11 ) 


where oti = gliOi)/{An) and 5'i(0) is the U(l) gauge coupling constant at the 
GUT scale. The corrections to the sparticle masses for the hrst two generations 
are listed in Appendix A. We note that although the non-universalities at the 
GUT scale in the hrst two generations vanish there are ehects at the electro- 
weak scale in the sparticle masses of the hrst two generations because of non¬ 
universalities in the Higgs sector and in the third generation via the trace 
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anomaly term. To get an idea of the size of the effects from this term, for 
Mg = GeV and ac = 1/24 one has p=0.0446. Thus the effect of the 
trace term is order a few percent for |5j| < 1. 

The analysis of non-universalities on the squark masses in the third gen¬ 
eration is more complicated. Here rri^ and toq obey the coupled renor¬ 
malization group equations (61)-(63) of Appendix A. Solution to these gives 
for rn\j^ 

3 

+ m\/2e{t) + Aomi/2f{t) -h mlh{t) - k{t)Al - -Sop (12) 
where Ah 2 is given by 


^\ s 2 + 6, + 6,) + 62 

and Dq dehnes the top Landau pole, i.e., 

yt n _ 1 « 

E{t)Do^ ° 

Here yo = h‘f/{ATTy where ht is the top Yukawa coupling, and 

16 3 13 

E{t) = (1 + /33(«))»»(1 + (32(())*ai + ft («))“■■ 


(13) 


(14) 


(15) 


In Eq. (15) jSi = Q;j(0)&i/47r, where 0:^(0) = ckg are the gauge coupling constant 
co-efficients at the GUT scale (cki = (5/3)Q;y), bi are the one loop beta function 
co-efficients dehned by (&i,& 2 ,& 3 ) = (33/5,1,-3), and E{t) = //U(t)(it. The 
non-universality effects in the stop matrix enter via ttIq and mp where the 
stop mass matrix is 


2 

mf 

tL 


—mt{At + pctnl3) 


mt{At + p,ctn/3)\ 


2 
m~ 

iR 


and are given by 


and 


1 2 

= rriQ + + {- — -siv?‘6w)M‘^cos2(3 


= mp -|- -|- {-)sin‘^9wMzCos2(3 

o 


(16) 


(17) 

(18) 
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Here and obey the set of coupled equations given in Appendix A. With 
the assumption of non-universalities of Eqs. (7) and ( 8 ) one has for 

2 1 

+ -ml + -{mlh - kAl + miA^f + m|e) 

+^(|/3 + / 2 -^/i)mi-^AoP (19) 

4:71 3 15 25 

where 

Aq = —^(52 -|- -|- 54 ) + Ss ( 20 ) 

Similarly for m|. one has 



+^{-f3 — f2 + -fi)ml/2 + -SoP (21) 

where 

Aq = —-—^(^2 + 53 + (^ 4 ) + <^4 (22) 

In the above, the non-universalities are all contained in the quantities Ah2, 
Aq, Afj and in the corrections involving the term S'op. We note that because 
of a peculiar accident that the sum of the corrections proportional to S'op in 
the subsectors involving H 2 ,Q and U vanishes, the trace anomaly receives no 
top quark Landau pole enhancement. A proof of this result is given in Ap¬ 
pendix A. An analysis similar to the above holds in the bottom squark sector. 
Details are given in Appendix A. 

In addition to above, the quantities p and are also effected. Using the 
radiative breaking relation Eq. (2) we hnd the following closed form solution 
for /i^ to one loop order; 

/i^ = mlCi + AIC 2 + m\C '4 + mi ^ 06*4 — -M| + —yS'op (23) 

Here 

Cl = ^(1 - + ^3 + S,)e) (24) 

C2 = -^'=1 C3 = -^(J - C, = -^/ (25) 
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where t = tanf3 and where the functions e,f,g,k are as dehned in Eqs. (3) and 
(4). Similarly, for m\ one has 


X 5 ^2 

m\ = ttiqDi + AqD 2 + m\D 3 + miAoD^ — 2 ^z + (26) 


where 


Dl - 7T5-7(1--) + ^2 _ I “ "2--(c>2 + C>3 + Oi)) (27) 


t2-l 


f2 + 1 f2 I X f2 I X 

Bs = - e). -04 = -^/ (28) 

In the above we have assumed the existence of non-universalities on phe¬ 
nomenological grounds. Theoretically there are several possible sources from 
where such non-universalities can arise. One source of non-universalities is 
the general Kahler potential. In general the Kahler potential can have gen¬ 
erational dependent couplings which lead to non-universalities in the visible 
sector after the breaking of supersymmetry in the hidden sector. However, 
even if the couplings in the Kahler potential were generation blind and the 
soft SUSY breaking masses for the scalars were universal at a scale higher 
than the GUT scale one will have non-universalities at the GUT scale arising 
from renormalization group running of the soft SUSY breaking parameters. 
Thus, for example, if one has universality of the soft parameters at the string 
scale, the running of these parameters from the string scale down to the GUT 
scale will generate non-universalities at the GUT scale. Model calculations 
indicate that such running typically generates |(5i 2 I < 0.5, and a similar level 
of non-universality in the relevant mass spectra at the electro-weak scale is ex¬ 
pected. However, enhancement of non-universalities can occur under special 
circumstances. To show this enhancement we consider the expression for 
Here the non-universalities are all contained in Ci. If the dominant universal 
terms cancel, then the non-universal effects become large. This situation can 
be easily seen to arise for large tan/?. Here for the case when mt ~ 167 GeV, the 
universal part cancels since Dq 1/3 and thus the non-universality effects get 
enhanced. A similar enhancement of non-universalities also occurs from the Aq 
dependence, although in a somewhat different manner. Here the enhancement 
occurs when the residue of the Landau pole in Aq vanishes. All three quan¬ 
tities, uiq and depend on Aq in the combination kA^ — fAQmi/ 2 - 

Using the renormalization group equation for Aq 

Ao = ^ - ^mi/2] HQ = tE - F 

JJq r 
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(29) 



where Ar = At — mi/ 2 (-f ^2 — H^/F) with H 2 dehned in Ref. (12), one hnds 
that 


kAl 


fAorn,,, = 1(1 - D„)A + Id „(1 - (30) 


2 F 

The Landau pole contribution vanishes when Ar = 0 which occurs when 


At ~ 0.61mg 0 - We give a graphical description of this phenomenon in Fig. 
1 where /i is plotted as a function of At. For the specihc set of parameters 
chosen in Fig. 1 the residue of the Landau pole vanishes when At/rriQ ~ 0.7. 
One hnds that Landau pole effects are suppressed close to this value of At, 
and there is considerable dispersion among the three curves corresponding to 
the sets of values = 0 = ^ 2 , = 1 = —^2 and = —1 = 82 . However, away 

from the region of At/rriQ ~ 0.7 one hnds that the Landau pole ehects begin 
to dominate and diminish the ehect of non-universalities. 

4. Constraints on Dark Matter 

We explore the parameter space of supergravity unihed models under the nat¬ 
uralness constraints which we assume to imply 


mo < 1 TeV, mg < ITeV, tan(3 < 25 


(31) 


The constraints of electro-weak breaking and the condition that there be no 
tachyons then imply the following range for T^ : 


At/mt) > —0.5 


(32) 


Constraints also arise from imposition of the condition that there be no nearby 
minima with lower energies that break color and charge conservation^^ (CCB 
constraints), and from experimental lower limits on the SUSY mass spectra 
given by CDF, DO and LEP experiments. 

In addition to the above there is an important constraint that arises from 
the decay 6 —>• s -t- 7 . This decay proceeds via loop corrections and is thus 
sensitive to physics beyond the Standard Model(SM). For the case of the SM, 
the loop contribution involves a W-t exchange, while for the supersymmetric 
case one has additional contributions arising from W — H~, W — t, and Z — b 
exchanges. Thus SUSY contributions are a priori as important as the SM 
contributions, and so one expects 6 —>• s -|- 7 decay to act as an important 
constraint on the parameter space of supergravity models. The CLEO value 
for the branching ratio for B ^ Xg + j is 


BR{B X, 7 ) = (2.32 ± 0.57 ± 0.35) x 10 


-4 


(33) 
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If one combines the errors in quadrature then BR(B—= (2.32 ± 0.67) x 
10“^. The SM prediction gives BR[B —>• = (3.28 ±0.33) x 10“^ for vnt = 

174 GeV|^, ^ 1^. In supersymmety, theoretical analyses when subject to 
the current experimental limits, impose serious constraints on the parameter 
space of the theory. The parameter space after the above constraints are im¬ 
posed must still be subject to the constraint of relic density of Eq. (1). The 
quantity fl^o is computed theoretically in each point in the supergravity 
parameter space using the relation]^, ^ 


^ 2.48 X 10 


-11 


Tv 


T. 


N. 


1/2 


2.73/ J{xf) 


(34) 


where T^ is the current background temperature, Tj is the freezeout temper¬ 
ature, {T^o/T^Y is the reheating factor, Nf is number of massless degrees of 
freedom at freezeout, Xf = kTf/m^^, and 

J (xf) = [ dx { av ) (x)GeV~^ (35) 

Jo 


Here < av > stands for the thermal everage, where a is the the neutralino an¬ 
nihilation cross-section and v is the neutralino relative velocity. In the analysis 
of the relic density we have used the accurate method for its computation p0|, 
1^, [3^ , ^ which carries out a correct thermal averaging over the Z and the 
Higgs poles which appear in the s channel in the neutralino annihilation. The 
accurate method for the computation of relic density is important as it has 
signihcant effects in the event rate analysis. The relic density computed in the 
above fashion is then subject to the constraints of Eq. (1) which further limits 
the parameter space of the theory. 

5. Detection of Neutralino Dark Matter We discuss now the possibilities 
for the detection of the neutralino dark matter. Several techniques, both 
direct and indirect, have been discussed over the years for its detection.The 
direct method includes (i) scattering by nuclei |^-|Q in terestial detectors. 


(ii) scattering by bound electrons[^, (iii) use of overheated microbubbles [^] , 
while the indirect methods include, (iv) annihilation in the center of sun and 
earth IEgI and 


V 


annihilation in the halo of galaxies. For our analysis 
here we limit ourselves to consideration of case (i). 

The total interaction in the scattering of neutralinos by quarks is given by 


Leff = (Xi7''75Xi)[97m(-4lTl ± AnPnjq] + {xiXi){QCmgq) (36) 
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where Pr^l = (1 ± 75)/2- The part with terms Al,Ar is the spin depen- 
dent(SD) part while the remainder is the spin independent (SI) part. The 
total event rate can thus be written as follows: 


R — [Rsi + Rsd] 


Pxi 

Vxi 

_0.3GeVcm~^. 

32t)km/s 


events 
kg da 


(37) 


where is the local mass density of xi incident on the detector, and is 
the incident Xi velocity. Rsi, is given by 


Rsi 


-^ 

[Mat + 


and the spin dependent rate is given by 


(38) 


Rsd = 




[Mjv + mxJ 


rA^ J(J + 1) I Asd 


(39) 


where J is the nuclear spin and A is dehned by < iV | X) S'i | iV >= A < 

I J \ N >. ^ 45 / and Asd are the corresponding amplitudes. We note that 
for large Mat, Rsi ~ Mat while Rsd ~ I/TTat- These results imply that the 
spin independent scattering becomes more dominant as the nucleus becomes 
heavier. Both Rsi and Rsd contain theoretical and experimental uncertain¬ 
ties. A major source of uncertainty in Rsi arises from the uncertainty in the 
matrix elements of the operator niqqq between nucleon states. There can be 
as much as 50% uncertainty in the strange quark contributions!^^ leading to 
a 30% uncertainty in Rsd For the case of the major source of ambiguity 
arises from the matrix elements of the axial current between nucleon states, 
i.e., of Aq dehned by < nlqx^XbQl'^ >=2s(jAg, where is the nucleon spin 
4-vector. There exist two sets of determinations of Aq, an old determination 
using the EMC data^^ and a new determination using the SMC datafpO 


For most nuclei the difference is not major and in our analysis here we use the 
values of Aq from the newer determination [Q. In addition both rates have 
uncertainties due to the nature of the nuclear form factors. Thus theoretical 
predictions are accurate to within perhaps a factor of two. 


6. Event Rates with Universal and Non-universal Soft Breaking 

We give now an analysis of the event rates and draw a comparison between the 
case where one uses universal boundary conditions at the GUT scale for the soft 
SUSY breaking parameters and the case when on has non-universal boundary 
conditions at the GUT scale. We begin with a discussion of event rates when 
one includes non-universalities in the Higgs sector, i.e., when (js = 0 = ^ 4 . 
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From Eqs. (23) and (24) one can see that a positive and a negative S 2 
make a positive contribution to while a negative and a positive 62 make 
a negative contribution to /i^. If one limits and 62 so that |(fj| < l(i=l,2), 
then the case = 1 = —62 gives the largest positive contribution to while 
the case (5i = —1 = —S 2 gives the largest negative contribution to /i^. These 
cases represent the extreme limits of how large the non-universality efffects 
can be within our prescribed limits on 5*. With the above in mind we focus 
on three cases to get an idea of the differences in the event rates between the 
universal and the non-universal cases: 

(i) (5i = 0 = (52, 

(ii) (5i = —1 = —S 2 , and 
(hi) (5i = 1 = -S 2 

In Fig. 2 we exhibit the maximum and minimum of event rates for xenon for 
the three cases listed above for fi > O.The analysis is done by mapping the 
full parameter space limited only by the naturalness constraints of Eq. (29), 
the constraints on the sparticle masses from CDF, DO and LEP experiments, 
the relic density constraint of Eq. (1) and the experimental constraint from 
& —> s -|- 7 of Eq. (31). The dips in the region below < 65 GeV arise 
due to the rapid annihilation in the vicinity of the Z pole and the Higgs pole. 
We note that the accurate method mentioned earlier for the computation of 
the relic density is important for the correct analysis of the event rates in this 
region. Comparison of cases(i) and (ii) shows that the event rates of (i) below 
the region < 65 GeV can be enhanced by a factor of 0(10) or more in 
the minimum curves for case(ii).This effect can be understood as follows: for 
case(ii) one has (5i < 0 and 82 > 0, which according to Eq. (23) drives 
towards the limit < 0 eliminating such points from the parameter space. 
Now this effect is more important for small tan/5. Since small tan/5 governs 
the minimum event rates one hnds that in this case since small tan/5 tend to 
get eliminated one has the effect that the lower limit of the event rates are 
raised. In contrast, for case (hi) one has (5i > 0 and (52 < 0 which gives the 
opposite effect as expected, i.e., on hnds that the minimum event rates are 
reduced by a factor of 0(10). However, above > 65 the non-universalities 
do not dramatically affect the event rates and the event rates have a uniform 
fall off as a function of the neutralino mass for all the three cases. A similar 
effect occurs for larger values of the top quark mass. However, here one hnds 
that a larger value of the top mass gives a larger Landau pole contribution 
which diminishes the relative contribution of the non-universal terms. 

We note that the spin independent interaction contributes signihcantly 
more than the spin dependent intereaction. This result holds not only for a 
heavy target such as xenon but also for retatively light targets such as CaF 2 . 
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In Tables 1, 2 and 3 we give a comparison of event rates for six additional 
target materials , i.e., He, CaF 2 , GaAs, Ge, Nal, and Pb, for the nniversal and 
non-nniversal soft SUSY breaking in the Higgs sector. (Recall, however, that 
there are signihcant uncertainties discussed above in the predictions.) Phenon- 
mena similar to the ones discussed for xenon also appear for these cases. Many 
of these materials listed in Tables 1, 2 and 3 are already being used as targets 
in the dark matter detectors. We note that for all the targets considered the 
event rates span several decades in magnitude within the allowed domain of 
the parameter space. Thus a full exploration of the parameter space certainly 
poses a formidable challenge to experimentalists. 

The analysis we have given so far is for the case p > 0. As discussed earlier 
for the case p < 0 a large part of the parameter space is eliminated because 
of the 6 —s + 7 constraint. Gonsequently the maximum event rates for the 
/i < 0 cases are much smaller than for the p > 0 case. A comparison of the 
allowed range of event rates for xenon for the cases p > 0 and /i < 0 is given 
in Table 4 (/x > 0) and Table 5 (/i < 0). One hnds that the event rates for 
the /r < 0 case are typically 0(10“^ — 10“^) smaller than for the /r > 0 case. 
Another important phenonmenon is that for both signs of p the parameter 
space corresponding to At/mo < —0.5 is eliminated. The reason that positive 
values of At are eliminated arises from the fact that the lightest stop mass 
turns tachyonic due to Landau pole effects, since the residue of the Landau 
pole A|j = {At — O.Olrrig)^ becomes large. 

In the analysis of this section thus far we assumed universality in the third 
generation sector and examined only the effects of non-universality in the Higgs 
sector. However, as discussed in Sec. 3 the non-universalities in the Higgs 
sector and in the third generation sector are highly coupled. Therefore, we 
turn now to a discussion of what we might expect, from the non-universality 
in the third generation sector. To get an idea of the nature of corrections 
one might expect, we display below the numerical sizes of the non-nniversal 
corrections • For values Mg = GeV, ac = 1/24 and 

rrit = 175 GeV one has Dq ~ 0.27 and from Eqs. (13), (20) and (22) one finds 


~ 0.64^2 — 0.36(53 + <^ 4 ) 

(40) 

Ai^ ~ 0.8853 - 0.12(52 + 54 ) 

(41) 

~ 0.7654 - 0.24(52 + 53 ) 

(42) 


From Eq. (40) we find that a positive S 2 in Ah 2 can be simulated by a negative 
value of S 3 or a negative value of ^ 4 , and a converse situation occurs for a 
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negative value of 62 - The effects of ^3 and ^4 on are more complicated 

since ^ 3,^4 enter with opposite signs. The stop masses do not enter 

directly in any signihcant manner in the amplitudes in the neutralino mass 
range we investigate, although they do affect the event rates by affecting the 
parameter space. In spite of this complexity we can glean the pattern in which 
( 53,(54 affect the event rates by examining the way they enter Ah 2 - In Fig. 3 
we exhibit the event rates for xenon when 61 = 62 = 64 = 0 and 63 takes 
on a range of values. The solid curve in Fig. 3 is for the case (53 = 0, the 
dotted curve for the case ^3 = 1, and the dashed for the case ^3 = — 1. Fig. 
3 shows that the minimum event rates in the region below < 65 GeV 
are lower relative to universal case for the case (53 > 0. This is as expected 
from our general analysis above, i.e., it is similar to the case ^2 < 0 in Fig. 2. 
Similarly the minimum event rates in the region < 65 GeV are enhanced 
relative to the universal case for ^3 < 0 which parallels the case ^2 > 0 for 
Fig. 2. This is again as expected. An analysis similar to the above for the 
case when (5i = (52 = (53 = 0 and 64 takes on a range of values is given in 
Fig. 4. The result of non-universalities is qualitatively similar to that of Fig. 
3 as expected. We note, however, that the parameter space allowed by the 
non-universalities in the third generation would be somewhat different and so 
one does not have a complete correspondence between the non-universalities 
in the Higgs sector (i.e., ^ 2 ) and those in the third generation sector (i.e., ^3 
and ^ 4 ). 

7. Effects of Grand Unification On Non-Universal Parameters 

In the previous section, the parameters describing non-universal soft SUSY 
breaking masses, i.e., 5 i,.., 54 , were viewed as independent quantities. This 
would be the case if the scale of SUSY breaking occured below Mg, since there 
only the SM gauge group is assumed to hold. However, if the SUSY breaking 
scale occurs above Mc^e.g., at the string scale Mstr), then the soft breaking 
masses at Mg must obey constraints imposed by the grand unihcation group 
G. We consider here this situation for the case where G contains SU(5) as a 
subgroup. This includes G=SU(5), SO(N) for N> 10 or E( 6 ) which essentially 
include all the grand unihcation groups that have been examined in signihcant 
detail. 

For this class of groups, we may characterize the light matter at Mg in 
terms of the SU(5) qunatum numbers, i.e., matter occurs in 10^ and 5^ repre¬ 
sentations (a=l,2,3 is a generation index) and the light Higgs doublet in a 5 
and 5 representation. One has then for the sfermionsf^ 


(?a — ^^La, diLa), Uq — Ca — ('Ra) 

5a (G = {yLa, ^La), da = djla') 


(43) 

(44) 
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The soft breaking masses have the form 


™10a = ™o(l + C); = ™0(1 + €) 


(45) 


with the Higgs masses obeying Eq. (7). There are therefore 8 independent soft 
breaking parameters for this general case. ( Eqs. (7, 43) contains 9 parameters 
6 ^, hi, 2 , ^ 0 , but actually one is redundant.). Thus we have the relations 

s,. = K = Se. = (46) 


For higher gauge groups, there can be additional relations. (For example, 
for SO(IO) with non-universalities in the third generation, some models of 
symmetry breaking give[^ h^=h 2 -hi.) 

One may impose further group constraints on the above parameters. Thus 
the suppression of FCNC suggests a near degeneracy in the hrst two gen- 
eraitons, which occurs naturally in a SU( 2 )i:i model where these generations 
are put into a doublet (d) representation of S\J{2)h and the third generation 
into a singlet (s) representation. Then Eqs. (45) reduce to 


= "io(l + <5(°)); mlo^=ml{l + 6l^)),a = l,2 (47) 

ml = ml{l + ml^ = ml{l + a = 1,2 (48) 


where quantities without generation indices are 811(2)/^ singlets. This model 
would have 6 independent parameters. 

While the FCNC constraints do not imply that the first two generation soft 
mass of the 10 representation is degenerate with that of the 5 representation, 
a further simplihcation occurs if one assumes that this is the case, i.e., one 
lets 

(49) 

One can then absorb the (1 -|- S(^d)) factor into the ml as the common mass for 
a=l,2, and use that as the reference mass for the third generation and Higgs 
masses, reducing the number of parameters to 5: , 5^^^, 5 i ,2 and mo. In the 

notation of Sec. 7, this is equivalent to non-universalities specihed by 


S 3 — 64 ] 5i; ^ 2 ; S 5 (50) 

where ^5 = = Si^ = Si,^. The parameter ^5 does not enter directly into 

the event rate R calculation or into the Landau pole or 6 —>■ s -f 7 analyses. 
It does, however, affect R indirectly in that it will enter into the relic density 
calculation and affect the determination of the parameter space that obeys 
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Eq. (1). However, there it enters only into the t-channel poles which generally 
make small contributions to the Xi annihilation cross section. We discuss now 
the effect of the model of Eq. (50) on event rates. For illustration we consider 
two cases with (5i(i=3,4) chosen at the extreme values under the reasonable 
constraint 6 i < 1. These are case(a): (53=(54=1, and case(b): (53=(54=-l(and (Js 
set to zero). For case(a) we hnd from Eqs. (13), (20) and (22) using Dq ~ 0.27 
that 

Ah^ ~ -0.73 + 0.6452 (51) 

Aq ~ 0.76 - 0.1252 (52) 

Ap ~ 0.51 - 0.2452 (53) 

For 52 in the range (1,-1) one hnds that Ah 2 lies in the range (-0.1,-1.37) while 
Aq lies in the range (0.64, 0.88) and Ap lies in the range (0.27,0.76). Because 
of the non-universalities arising from 53 = 54 > 0 this case is again different 
from the cases we have discussed before. However, we can still draw some 
rough comparisons with the previous cases. Thus, for example, we note that 
since the effective value of 62 in this case is negative the minimum of the event 
rates should be similar to the minimum dotted curve of Fig. 2 where 52 was 
also negative. The result of the analysis is presented in Fig. 5. We see that 
indeed the minimum of the event rates here is similar to the minimum of the 
event rates in Fig. 2 as expected. In general, the effect of non-universalities 
is not very large here. A similar analysis holds for case(b). Here again using 


Eqs. (13), (20), (22) and Dq ~ 0.27 we hnd that 

Ah 2 ^ 0.73 + 0.6452 

(54) 

Aq ~ -0.76 - 0.1252 

(55) 

Ap ~ -0.51 - 0.2452 

(56) 


For 52 in the range (1,-1) one hnds Ah^ in the range (0.1,1.37), Aq in the 
range (-0.64, -0.88), and Ap in the range (-0.27,-0.76). We note that these 
ranges are exact mirrors of the ranges for case (a) but with a negative sign. 
Of course, the full dynamics of this system is very complicated since the al¬ 
lowed parameter space is affected differently from previous cases because of 
the specihc nature of non-universalities here. However, as in case(a) we can 
glean some understanding by a comparison with the previous analyses. We 
note that the effective value of S 2 is positive here. Thus one expects that the 
minimum event rates might show the same pattern as the positive 52 case of 
Fig. 2. The result of the analysis is exhibited in Fig. 6. We hnd that indeed 
the minima curves of Fig. 6 have a behavior similar to the dashed curve of 
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Fig. 2 as expected. The remarkable feature here is the increase by a factor 
of 10 in the maximum event rates when S 2 is positive. This is due to the fact 
that 62 , 63 , 64 ^ now act together to reduce (see Eqs. (23) and (24)) which 
increases the maximum event rates^®. 

8. Effects of More Accurate Determination 

As mentioned in the Introduction the next round of satellite experiments are 
expected to determine the Hubble parameter to within 0(10%) or better accu¬ 
racy. Such an accurate determination of the Hubble parameter will put more 
stringent bounds on the neutralino relic density and hence more strin¬ 

gently constrain the SUGRA parameter space. We investigate here the effects 
of these more stringent constraints on the neutralino dark matter event rates. 
For specihcity we consider a narrow band of in the range 

0.225 < < 0.275 (57) 

Fig. 7 gives the maximum and the minimum event rates as a function of the 
Xi mass. A comparison of Figs. 2 and 7 shows that a narrower range 
more sharply constrains the SUGRA parameter space and the allowed band 
of event rates shrinks. The allowed range of the neutralino mass is also more 
sharply constrained since the mass is constrained to lie below ~ 90 GeV, 
as above this mass exceeds the upper limit of Eq. (57). Since the neu¬ 
tralino mass and the gluino mass scale one also deduces an upper limit on 
the gluino mass of about 650 GeV. Thus a precise measurement of will 

put an interesting upper bound on the gluino mass which would be testable 
at accelerators such as the upgraded Tevatron and the LHG. In general the 
overall effect of the more stringent constraints on h? is to eliminate a part 
of the parameter space allowed by Eq. (1) and thus limit more severely the 
permissible domain of the event rates. 

9. Conclusions 

In this paper we have investigated the effects of non-universalities in the Higgs 
secotor and in the third generation squark sector on the low energy param¬ 
eters. We have shown that the Higgs sector non-universalities and the non¬ 
universalities in the third generation squark sector are highly coupled and we 
have exhibited analytic forms of their effects on the parameters p, vriA and the 
third generation squark masses at low energy. Next we investigated the effects 
of these non-universalities on the event rates in the scattering of neutralinos 
by nuclei. Here we found that the minimum event rates in the neutralino mass 
region 30 GeV < < 65GeU can be increased or decreased by a factor 

of 0(10) depending on the sign of the non-universality. In contrast, in the 
neutralino mass region > 65 GeV the minimum and the maximum event 
rates are not appeciably affected. These results have important implications 
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for dark matter detectors. Supergravity unified models with universal bound¬ 
ary conditions for soft SUSY breaking parameters give event rates which lie in 
a wide range R= 0(1 — 10“®) events/kg da Thus while the current generation 
of dark matter detectors can sample a part of the parameter space, one needs 
far more sensitive detectors (more sensitive by a factor of 10 ^) to sample the 
entire parameter space of the model. However, inclusion of non-universalities 
shows that for certain signs of the non-universality, the minimum event rates 
can increase by a factor of 10-100 in the region < 65GeV. These models 
can be completely tested in this mass region if the detector sensitivity could 
be enhanced by a factor of about 10^. Such an enhancement does not seem 
out of reach with new techniques currently being discussed [^]. Another inter¬ 
esting result seen was the appearance of large event rates, up to 10 event/kg 
da for the case when ^3 = ^4 = —1. This case looks very interesting from 
the view point of experimental detection of dark matter. We have also anal¬ 
ysed the implications of more stringent constraints on Qh^ which is expected 
from the next round of satellite experiments, and found that while the event 
rates typically become smaller a signihcant part of the parameter space yields 
event rates which are within reach of dark matter experiments currently being 
planned. 
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Appendix A 

All the scalar sparticle masses receive contributions from the Tr(Ym^) term. 
Corrections to the first two generations of sparticle masses are given by 

( 58 ) 

= -^Sop (59) 

where So and p are dehned in the text. Under the assumption that the squarks 
and slepton masses are universal at the GUT scale for the hrst two generations, 
but with non-universaliites in the Higgs sector and in the third generation we 
can deduce the following sum rule for the Tr(Ym^) term: 

pTr{Yw?)„ = -^sirfewCOs2l3Ml + 

( 60 ) 

Here Tr{Ym?)o is the value of TriYmS) at the GUT scale while the right 
hand side is evaluated at the electro-weak scale. We see from above that if 
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the masses of the first two generations of squarks and sleptons are known with 
enough precision, one can determine the value of the trace anomly term at the 
GUT scale. For the case when Tr{Yw?)o vanishes the above sum rule reduces 
to the one given in Ref. (56). 

satisfy the following set of coupled equations 

d/TYt^ 3 

—7^ = “SFtSt + (3Q;2iU2 + —diiVf^) + 'yiOiiS (61) 

at 5 

dm% T 6 o 16 Ox , X 

+ {-a,Mi + -a,M^) + ^ 2^18 (62) 

at 6 15 

dm% 16 „ o 1 „ 

—— = —YtTif + + 3(52^2 + (63) 

at 3 15 

where + A^), and 7i=-^, 72 =| and 73 =-^ and S is 

the running parameter with S(0)=S'o. We show now that the evolution of the 

trace anomaly term is not affected by the top Yukawa coupling. To see this it 

is best to go to a decoupled set of equations using the functions ei(t) given by 


^'h 2 ~ + 62+63 

(64) 

= 2 ei — 2 e 2 

(65) 

mQ = 61 + 62 - 63 

( 66 ) 


It is easily seen that the evolution of 62 , and 63 does not involve the Yukawa 
coupling and only the evolution of ei does. One has 

^ = -6Y*ei-Y*R? + (^a3M32+a2M2 + gaiM2) + l(7i+72+73)di^ (67) 
However, since 

(71 + 72 + 73) = 0 (68) 

one hnds that the evolution of the anomaly terms in Eqs. (61)-(63) are not 
affected by the top Yukawa coupling and are thus not sensitive to the Landau 
pole effects. The bottom squark mass matrix is given by 


where 



(69) 

(70) 
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and where is defined in the text. 

Figure Captions 

Fig.l: /i vs At for the case when = 0 = ^4, mo = 300 GeV, rrig = 350 
GeV, and (i)5i = 0 = ^2(solid), (ii) (5i = —1 = —52(dashed), and (iii)5i = 1 = 
—^2(dotted). Here vanishes for At/rriQ ~ 0.7. 

Fig. 2: Maximum and minimum of event rates/kg da for xenon for /i > 0 for 
the case when ^3 = 0 = ^4 and (a)5i = 0 = ^2(solid), (b)5i = 1 = —^2(dotted), 
and (c)5i = —1 = —82 (dashed) when 0.1 < < 0.4, and mi=175 GeV. 

Fig. 3: Maximum and minimum of event rates/kg da for xenon for the case 
when 5i = ^2 = ^4=0 and (a)53 = O(sohd), (b)53=l(dotted), and (c)53=-l 
(dashed) when 0.1 < < 0.4, and mt=175 GeV. 

Fig. 4: Maximum and minimum of event rates/kg da for xenon for the case 
when 5i = ^2 = ^3=0 and (a)54 = O(sohd), (b)54=l(dotted), and (c)54=-l 
(dashed) when 0.1 < < 0.4, and mt=175 GeV. 

Fig. 5: Maximum and minimum event rates when ^3 = ^4 = 1, ^5 = 0 and 
(a) 5i = 0 = 52(solid), (b) 5i=l= -52 (dotted), and 5i=-l=-52 (dashed) when 
0.1 < < 0.4, and mt=175 GeV. 

Fig. 6: Same as Fig. 5 when 53 = 54 = —1. 

Fig. 7: Maximum and minimum of event rates/kg da for xenon for the case 
when 53 = 54=0 and (i) 5i = 0 = 52(sohd), (ii)5i = —1 = —52(dashed), and 
(iii)5i = 1 = —52(dotted) when 0.225 < < 0.275, and mi=175 GeV. 
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41.5 

41.5 

52.4 

52.4 


RmIti 

RMax 

RmItl 

Rlilax 

He 

6.97E-4 

3.43E-3 

6.89E-4 

2.49E-3 

GaF2 

8.33E-4 

4.26E-3 

9.02E-4 

2.25E-2 

GaAs 

1.46E-4 

1.40E-3 

1.58E-4 

8. MEM 

Ge 

1.08E-4 

1.23E-3 

l.llE-4 

8.33E-2 

Nal 

1.18E-4 

1.68E-3 

1.14E-4 

0.12 

Pb 

1.82E-4 

2.83E-3 

1.68E-4 

0.21 

(GeV) 

61.4 

61.4 

71.2 

71.2 


RmIu 

Rhlax 

RmIti 

RjHax 

He 

2.95E-5 

9.35E-3 

3.05E-7 

4.28E-3 

GaF2 

3.53E-4 

5.48E-2 

1.51E-4 

2.45E-2 

GaAs 

1.32E-3 

0.17 

6.41E-4 

7.89E-2 

Ge 

1.35E-3 

0.18 

6.57E-4 

8.04E-2 

Nal 

2.02E-3 

0.27 

9.65E-4 

0.11 

Pb 

3.38E-3 

0.45 

1.58E-3 

0.19 

(GeV) 

80.7 

80.7 

117.5 

117.5 


f^Min 

Rlvlax 

Rmiti 

Rlvlax 

He 

5.23E-6 

2.31E-3 

l.OlE-4 

1.87EM 

GaF2 

8.95E-5 

1.27E-2 

2.29E-4 

8.19EM 

GaAs 

3.47E-4 

3.99E-2 

3.01E-4 

2.22E-3 

Ge 

3.55E-4 

4.06E-2 

2.96E-4 

2.25E-3 

Nal 

5.MEM 

5.86E-2 

4.03E-4 

3.11E-3 

Pb 

8.32E-4 

9.45E-2 

6.27E-4 

4.86E-3 


Table 1: Minimum and maximum of event rates/kg da are given for the case 
5i = 0 = (52 and /i > 0 for several neutralino masses for the target materials 
He, CaF2, GaAs, Ge, Nal, and Pb. 
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(GeV) 

41.6 

41.6 

53.9 

53.9 


RlVlin 

RjVlax 

RjVlin 

RjVlax 

He 

3.57E-5 

3.6IE-5 

7.88E-6 

1.91E-3 

CaF2 

4.48E-5 

1.034E-3 

1.07E-5 

1.79E-2 

GaAs 

1.62E-5 

4.10E-3 

3.05E-6 

6.53E-2 

Ge 

1.45E-5 

4.21E-3 

2.53E-6 

6.68E-2 

Nal 

1.99E-5 

6.56E-3 

3.20E-6 

0.10 

Pb 

3.38E-5 

1.15E-2 

5.13E-6 

0.17 

m^^{GeV) 

61.0 

61.0 

70.1 

70.1 


RlS/lin 

^Max 

RlS/lin 

RjS/lax 

He 

3.18E-4 

1.09E-2 

1.48E-4 

5.05E-3 

GaF2 

6.96E-4 

6.06E-2 

3.52E-4 

2.77E-2 

GaAs 

1.16E-3 

0.19 

6.30E-4 

8.81E-2 

Ge 

1.17E-3 

0.19 

6.32E-4 

8.98E-2 

Nal 

1.72E-3 

0.29 

9.18E-4 

0.13 

Pb 

2.86E-3 

0.49 

1.50E-3 

0.21 

m^^{GeV) 

80.9 

80.9 

117.0 

117.0 


f^Min 

RjVlax 

^Min 

^Max 

He 

1.02E-5 

2.15E-3 

l.lOE-4 

1.80E-4 

GaF2 

8.87E-5 

1.14E-2 

2.28E-4 

7.85E-4 

GaAs 

3.14E-4 

3.53E-2 

2.36E-4 

2.12E-3 

Ge 

3.20E-4 

3.60E-2 

2.29E-4 

2.14E-3 

Nal 

4.63E-4 

5.19E-2 

3.09E-4 

2.97E-3 

Pb 

7.49E-4 

8.37E-2 

4.78E-4 

4.63E-3 


Table 2: Same as Table 1 except that (5i = 1 = —82 
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(GeV) 

50.3 

50.3 

54.6 

54.6 


RmItl 

^Max 

RmItl 

^Max 

He 

6.61E-3 

1.50E-2 

1.14E-4 

4.36E-2 

CaF2 

8.90E-3 

0.10 

7.33E-4 

0.28 

GaAs 

3.11E-3 

0.36 

2.47E-3 

0.95 

Ge 

2.72E-3 

0.37 

2.52E-3 

0.97 

Nal 

3.64E-3 

0.57 

3.82E-3 

1.47 

Pb 

5.99E-3 

0.97 

6.47E-3 

2.49 

(GeV) 

61.2 

61.2 

70.9 

70.9 


RmIti 

Rhlax 

RmIti 

Rhlax 

He 

3.71E-5 

2.80E-2 

3.74E-6 

7.30E-2 

GaF2 

4.09E-4 

0.14 

1.81E-4 

0.16 

GaAs 

1.52E-3 

0.45 

7.46E-4 

0.29 

Ge 

1.55E-3 

0.46 

7.65E-4 

0.29 

Nal 

2.32E-3 

0.68 

1.12E-3 

0.43 

Pb 

3.88E-3 

1.14 

1.84E-3 

0.71 

m^^{GeV) 

80.5 

80.5 

117.5 

117.5 


^Min 

^Max 

^Min 

F^Max 

He 

1.65E-6 

8.20E-2 

1.04E-4 

8.02E-2 

GaF2 

9.40E-5 

0.18 

2.37E-4 

l.lOE-2 

GaAs 

3.88E-4 

0.33 

3.12E-4 

1.89E-2 

Ge 

3.98E-4 

0.33 

3.07E-4 

1.23E-2 

Nal 

5.76E-4 

0.49 

4.18E-4 

1.14E-2 

Pb 

9.33E-4 

0.81 

6.50E-4 

1.57E-2 


Table 3: Same as Table 1 except that (5i = — 1 = —82 
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At 

5i — 0 — ^2 

5i — 1 — —^2 

5i — —1 — —62 

-0.5 

Rmiti Rhlax 

RMin Rhdax 

3.88E-6 2.83E-3 

RmIu Rhiax 

0 

5.97E-5 1.28E-2 

6.43E-6 1.53E-2 

4.40E-3 0.19 

1.0 

2.28E-3 0.26 

l.lOE-5 7.14E-2 

1.15E-2 1.82 

2.0 

2.73E-3 0.36 

1.80E-3 0.20 

5.14E-3 1.22 

3.0 

6.83E-4 0.52 

5.82E-4 0.36 

8.14E-4 0.79 

4.0 

3.19E-4 0.53 

2.90E-4 0.43 

3.52E-4 0.84 

5.0 

3.10E-4 0.56 

3.53E-4 0.45 

3.42E-4 0.71 

6.0 

2.14E-4 0.46 

3.84E-4 0.37 

2.22E-4 0.58 

7.0 

1.82E-4 0.21 

1.66E-4 0.17 

2.01E-4 0.25 


Table 4: A comparison of the minimum and maximum event rates/kg da for 
xenon for the cases (5i = 0 = 52, 5i = 1 = — 62 , and 5i = —1 = — 62 , and /i > 0 
as a function of Af. 


At 

5i — 0 — ^2 

5i — 1 — —82 

5i — —1 — —82 

-0.5 

0.0 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

R-Min Rjvlax 

4.65E-6 6.40E-3 
4.65E-6 8.15E-3 
4.66E-6 1.48E-3 
4.67E-6 6.18E-4 

4.67E-6 5.55E-4 
4.67E-6 4.40E-4 
4.67E-6 3.49E-4 
4.67E-6 3.09E-4 
4.67E-6 3.17E-4 

RmIu Rjvlax 

2.55E-6 1.33E-2 
2.22E-6 2.33E-2 
4.67E-6 2.65E-3 
4.67E-6 3.65E-3 
4.67E-6 5.49E-4 
4.67E-6 4.36E-4 
4.67E-6 3.46E-4 
4.67E-6 2.77E-4 
4.67E-6 3.15E-4 

Rmiu Rhdax 
4.65E-6 1.24E-4 
4.65E-6 1.73E-2 
4.67E-6 9.04E-4 
4.67E-6 6.21E-4 

4.67E-6 5.60E-4 
4.67E-6 4.44E-4 
4.67E-6 l.OlE-3 
4.67E-6 3.11E-4 
4.67E-6 l.OE-2 


Table 5: Same as Table 4 except that /i < 0. 


24 





References 

[1] P.F. Smith and J.D. Lewin, Phys.Rep.187,No.5,203(1990). 

[2] See,e.g., J. Primack et al. Phys. Rev. Lett. 74, 2160 (1995). 

[3] G. Jungman, M. Kamionkowski and K. Greist, Phys. Rep. 267,195(1995). 

[4] V. Berezinsky, A. Bottino, J. Ellis, N. Forrengo, G. Mignola, and S. 
Scopel, Astropart. Phys.5:1-26(1996). 

[5] A.H. Ghamseddine, R. Arnowitt and P. Nath, Phys. Rev. Lett. 29, 970 
(1982). For reviews see P.Nath, R. Arnowitt and A.H.Ghamseddine, “Ap¬ 
plied N = 1 Supergravity” (World Scientihc, Singapore, 1984); H.P. Nilles, 
Phys. Rep. 110, 1 (1984); R. Arnowitt and P. Nath, Proc. of VII J.A. 
Swieca Summer School ed. E. Eboli (World Scientihc, Singapore, 1994). 

[6] R. Barbieri, S. Ferrara, and Savoy, Phys. Lett. B119,343(1982). 

[7] L.Hall, J.Lykken and S.Weinberg, Phys. Rev. D27,2359(1983). 

[8] P. Nath, R. Arnowitt and A.H. Ghamseddine, Nucl. Phys. 
B227,121(1983). 

[9] K. Inoue et al. Prog. Theor. Phys 68, 927 (1982); L. Ibanez and G.G. 
Ross, Phys. Lett. BllO, 227 (1982); L. Avarez-Gaume, J. Polchinski 
and M.B. Wise, Nucl. Phys. B221, 495 (1983), J. Ellis, J. Hagelin, D.V. 
Nanopoulos and K. Tamvakis, Phys. Lett. B125, 2275 (1983); L.E. Ibanez 
and G. Lopez, , Nucl. Phys. B233, 545 (1984); 

[10] J. Ellis, G. Ridolh and F. Zwirner, Phys. Lett. B262, 477 (1991). 

[11] R. Arnowitt and P. Nath,Phys.Rev.D46,3981 (1992). 

[12] L. Ibanez, G. Lopez, and G. Munos, Nucl. Phys.B256, 218(1985). 

[13] S.P. Martin and M.T. Vaughn, Phys. Rev. D50, 2282(1994). 

[14] D. Pierce and A. Papadopoulos, Nucl. Phys. B430,278(1994). 

[15] R. Arnowitt and P. Nath, Phys.Rev. Lett. 69,725(1992); P. Nath and R. 
Arnowitt,Phys.Lett. B289,368 (1992); 

[16] R. Arnowitt and P. Nath, Phys. Rev. D54,2374(1996). 


25 



[17] S. K. Soni and H. A. Weldon, Phys. Lett. B126, 215(1983);V. S. 
Kaplnnovsky and J. Lonis, Phys. Lett. B306, 268(1993). 


[18] D. Matalliotakis and H.P. Nilles, Nncl.Phys.B435, 115(1995). 

[19] M. Olechowski and S. Pokorski, Phys.Lett. B344, 201(1995). 

[20] N. Polonski and A. Pomarol,Phys.Rev.D51,6532(1995). 

[21] P. Nath, J. Wn and R. Arnowitt, Phys.Rev.D52,4169(1995). 


[22] J.M. Frere, D.R.T. Jones and S.Raby, Nncl. Phys.B222, 11(1983); 
M.Drees, M.Glnck and K.Grassie, Phys. Lett.B157, 164(1985); J.A. 
Gasas, A. Lleyda and G. Mnnoz, Nncl. Phys. B471, 3(1996); H.Baer, 
M. Brhlik. D. Gastano, hep-ph/960801; see however, A. Knsenko, P. Lan- 
gacker, and G. Segre, Phys. Rev.D54, 5824(1996); A. Riotto, E.Ronlet 
Vilja, |hep-ph/96074031 


and L. 
ph/960834q. 


A.Knsenko, and P. Langacker, |hep- 


[23] M.S. Alamet ah (GLEO Gollaboration), Phys. Rev. Lett. 74, 2885 (1995). 


[24] S. Bertolini, F. Borznmati and A. Masiero, Phys. Rev. Lett. 59, 180 
(1987); B. Grinstein, R. Springer and M.B. Wise, Nncl. Phys. B339, 
269 (1990); S. Bertolini, F. Borznmati, A. Masiero and G. Ridolh, Nncl. 
Phys. B353, 591 (1991); R. Barbieri and G. Gindice, Phys. Lett. B309, 
86 (1993); M. Misiak, Phys. Lett. B269, 161 (1991); Nncl. Phys. B393, 23 
(1993); F. Borznmati, M. Drees and M. Nojiri, Phys.Rev. D51,341(1995). 


[25] A.J. Bnras, M. Misiak, M. Miinz and S. Pokorski, Nncl. Phys. B424, 374 
(1994); M. Ginchini, E. Franco, G. Martinelli, L. Reina and L. Silvestrini, 
Phys. Lett. B316, 127 (1993). 


[26] M. Misiak, |hep-ph/9612313| . 


[27] P. Nath and R. Arnowitt, Phys.Lett.B336, 395(1994);H. Baer and M. 
Brhlick, |hep-ph/9610224| . 


[28] B.W. Lee and S. Weinberg, Phys. Rev. Lett. 39, 165 (1977); D.A. Dicns, 
E. Kolb and V. Teplitz, Phys. Rev. Lett. 50, 168 (1977); H. Goldberg, 
Phys. Rev. Lett. 50, 1419 (1983); J. Ellis, J.S. Hagelin, D.V. Nanoponlos 
and N. Srednicki, Nncl. Phys. B238, 453 (1984); J. Lopez, D.V. Nanopon¬ 
los and K. Ynan, Nncl.Phys.B370,445(1992); M. Drees and M.M. Nojiri, 
Phys. Rev. 047,376(1993). 


26 












[29] For a review see E.W. Kolb and M.S. Turner, “The Early Universe” 
(Addison-Wesley, Redwood City, 1989). 

[30] K. Greist and D. Seckel, Phys. Rev. D43, 3191 (1991); P. Gondolo and 
G. Gelmini, Nucl. Phys. B360, 145 (1991). 

[31] R. Arnowitt and P. Nath, Phys. Lett. B299, 58 (1993); B303, 403 (1993) 
(E). 

[32] P. Nath and R. Arnowitt, Phys. Rev. Lett. 70, 3696 (1993). 

[33] H. Baer and M. Brhlick, Phys.Rev.D53,597(1996). 

[34] M.W. Goodman and E. Witten, Phys. Rev. D31, 3059(1983); K. Greist, 
Phys. Rev. D38, (1988)2357; D39,3802(1989)(E); J. Ellis and R. Flo¬ 
res, Phys. Lett. B300,175(1993); R. Barbieri, M. Frigeni and G.F. Gin- 
dice, Nucl. Phys. B313,725(1989); M. Srednicki and R.Watkins, Phys. 
Lett. B225,140(1989); R. Flores, K. Olive and M. Srednicki, Phys. Lett. 
B237,72(1990). 

[35] A. Bottino etal. Astro. Part. Phys. 1, 61 (1992); 2, 77 (1994). 

[36] V.A. Bednyakov, H. V. Klapdor-Kleingrothaus and S.G. Kovalenko, Phys. 
Rev. 050,7128(1994). 

[37] R. Arnowitt and P. Nath, Mod. Phys. Lett. A 10,1257(1995). 

[38] P. Nath and R. Arnowitt, Phys. Rev. Lett.74,4592(1995); R. Arnowitt 
and P. Nath, Phys. Rev. 054,2374(1996). 

[39] E. Diehl, G. Kane, G. Kolda and J. Wells, Phys. Rev.052,4223(1995). 

[40] L. Bergstrom and P. Gondolo, |hep-ph/9510252| . 

[41] J.D. Vergados, J. Phys.G22,253(1996). 

[42] G.L. Kane and J.D. Wells, Phys.Rev. Lett.76, 4458(1996). 

[43] B. de Garlos and G.V. Kraniotis, Piep-ph/9610355| . 

[44] G.D. Starkman and D.N. Spergel, Phys.Lett.74,2623(1995). 

[45] J.I. Gollar, Phys.Rev.D54,1247(1996). 

[46] M. Kamionkowski, Phys.Rev.O44,3021(1991); F. Halzen, T. Saltzer and 
M. Kamionkowsi, Phys.Rev. 045,139(1992). 


27 






[47] M. Kamionskowski, K. Greist, G. Jungman, and B. Sadoulet, 
Phys.Rev.Lett. 74,5174(1995). 

[48] T.P. Gheng, Phys. Rev. D38, 2869 (1989); J. Gasser, H. Leutwyler and 
M.E. Saino, Phys. Lett. B253, 252 (1991). 

[49] J. Ashman etal (EMG Gollaboration), Nucl. Phys. B328, 1 (1989); A. 
Manohar and R. Jaffe, Nucl. Phys. B337, 509 (1990). 

[50] D. Adams et ah (SMG Gollaboration), Phys. Lett. B329, 399 (1994); 
R. Arnold et al. (E143 Gollaboration) in “Intersections of Particle and 
Nuclear Physics”, St. Petersburg, EL (1994); J. Ellis and M. Karliner, 
Phys.Lett.B341,397(1995). 

[51] For the flipped SU(5) model|^ one makes the changes u <->• d, e < -> 9 

with en appearing in an SU(5) singlet. 

[52] 1. Antoniadis, J. Ellis, J. Hagelin and D.V. Nanopoulos, Phys. Lett. B208, 
209(1988). 

[53] Y. Kawamura, H. Murayama and M. Yamauguchi, Phys. Lett. B304, 
52(1994); R. Arnowitt and P. Nath, [hep-ex/9605011| ;GTP-TAMU-02/97/ 
NUB-TH- 3552/97. 

[54] M. Dine, R. Leigh and A. Kagan, Phys. Rev. D48, 4269(1993). 

[55] D. Gline.”On a Discriminatory Liquid Xenon SUSY Wimp Detector”, 
Nucl. Phys. B (in press); P. Benetti et al, Nucl. Inst, and Method for 
Particle Physics Research, A307,203 (1993). 

[56] S. Martin and P. Ramond, Phys. Rev. D48, 5365(1993). 


















